Effect of reheating on electroweak baryogenesis 
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The latent heat released during the expansion of bubbles in the electroweak phase transition reheats 
the plasma and causes the bubble growth to slow down. This decrease of the bubble wall velocity 
affects the result of electroweak baryogenesis. Since the efficiency of baryogenesis peaks for a wall 
velocity ~ fO~^, the resulting baryon asymmetry can either be enhanced or suppressed, depending 
on the initial value of the wall velocity. We calculate the evolution of the phase transition taking 
into account the release of latent heat. We find that, although in the SM the baryon production is 
enhanced by this effect, in the MSSM it causes a suppression to the final baryon asymmetry. 
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Generating the baryon asymmetry of the universe 
(BAU) at the electroweak phase transition is a very at- 
tractive idea Q|. Although in the minimal Standard 
Model (SM) it fails to explain the observed baryon abun- 
dance, such a proposal can be successful in an exten- 
sion of the model. It has been shown, for instance, that 
electroweak baryogenesis is quantitatively possible in the 
Minimal Supersymmetric Standard Model (MSSM), pro- 
vided that the Higgs boson and the lightest stop are suf- 
ficiently light |,|. 

According to the standard mechanism baryogene- 
sis occurs near the walls of expanding bubbles that form 
during the phase transition. After a bubble is nucleated 
and begins to grow, its wall quickly reaches a terminal ve- 
locity in the hot plasma Due to CP violating inter- 
action of particles with bubble walls, different densities 
of left-handed quarks and their antiparticles are built up 
in front of the walls. This left-handed asymmetry biases 
the anomalous baryon number violating sphaleron inter- 
actions in the symmetric phase. As a consequence, a 
net baryon asymmetry is generated in front of the walls 
and immediately caught by them. In the broken sym- 
metry region inside the bubbles sphaleron processes are 
suppressed, so a subsequent washout of the baryon asym- 
metry is avoided. 

The generated BAU has a strong dependence on the 
bubble wall velocity. If the latter is too large, the left- 
handed density perturbation will pass so quickly through 
a given point in space that sphaleron processes will not 
have enough time to produce baryons. Thus the result- 
ing BAU will be small. On the other hand, for very 
small velocities thermal equilibrium will be restored, so 
the baryon asymmetry will be erased by sphalerons and 
the BAU will be small again. Consequently, the baryon 
production has a maximum at an intermediate wall ve- 
locity. Comparison of the baryon number violation time 
scale with the time of passage of the chiral asymmetry 
|^,||^ gives a wall velocity 10~^ for maximum 

baryon asymmetry. Recent numerical calculations con- 
firm that the BAU tends to peak for such a small value 



On the other hand, recent calculations of the friction 
of the plasma indicate that the wall velocity can be of 
that order of magnitude jljjl^. However, these esti- 
mates of Vw do not take into account the hydrodynam- 
ics of the phase transition ||9|,|T^. For such small veloc- 
ities the only effect of hydrodynamics is a homogeneous 
reheating of the plasma during the expansion of bub- 
bles. As is much less than the speed of sound in 
the plasma, the latent heat liberated in the expansion 
of a bubble is quickly distributed throughout space by 
a shock front that precedes the propagating wall of the 
bubble. As a consequence of this uniform reheating the 
bubble expansion slows down. It was shown by Heck- 
ler that the decrease of Vw can dramatically affect 
the result of electroweak baryogenesis. The importance 
of reheating can be estimated by comparing the la- 
tent heat L = T{d (AV) /dT)^^^,^ , where (T) is the 
free energy difference between the symmetric and bro- 
ken symmetry phases, with the energy needed to bring 
the plasma back to the critical temperature Tc from the 
temperature r„ at which nucleation of bubbles begins, 
Ap = (TT^g^/SO) (T^ - T^), where ~ 107 is the num- 
ber of degrees of freedom of the plasma. In the SM, L is 
at least one order of magnitude less than Ap, but in the 
MSSM the two quantities are of the same order |]l^ , |l5| , 
so the effect of reheating can be important in that case. 

In this paper we compute the evolution of the phase 
transition including this effect in a simple model with 
a one-Higgs effective potential. Using an analytical ap- 
proximation for the dependence of the BAU on Vw^ we 
evaluate the effect of the decrease of the bubble wall ve- 
locity on baryogenesis. It is known that in the SM the 
reheating enhances electroweak baryogenesis . On the 
contrary, since in the MSSM the wall velocity is initially 
in the range of maximum BAU, its decrease will cause a 
suppression to the baryon asymmetry. 

We use an effective potential of the form iTlplJT^ 



V ((/), T) = D (T^ - T^] 



b^-ETct>' + \<^\ (1) 
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which possesses a first-order phase transition between 
the critical temperature Tc = Tq/ yjl — E'^ /\D and Tq. 
When the universe cools below Tc ~ lOOGey, bubbles 
of the broken symmetry phase begin to nucleate with a 
probability per unit volume and time T ~ 7"4g-S3(T)/T^ 
where [4> {r)] is a three-dimensional instanton action 
that coincides with the energy of the nucleated bubble 
| p6[ . The configuration (r) of the bubble is obtained by 
finding an extremum of S'3 . It can be calculated numeri- 
cally to obtain the nucleation rate T and the radius tq 
and wall width 1^ of the nucleated bubble ||l5| as func- 
tions of temperature. 

For the terminal velocity of the bubble wall in the 
plasma we use the formula (see also |]5|,|6|,p|Jl5[| ) 



20r/^ (T) AV (T) 



(2) 



where 77 is a dimensionless friction coefficient accounting 
for the viscosity of the plasma, and v (T) is the minimum 
of V {4>, T) . The main effect of the reheating will be to 
decrease the free energy difference (T) and hence 
the value of Vw ■ The progress of the phase transition is 
determined by the fraction of space that is still in the 
symmetric phase |18|, 



/ {t) = exp 



47r 



r (T') r {t', ty dt' 



(3) 



where r {f, t) = tq (T') + J* {T") dt" is the radius at 
time i of a bubble created at t' . The variation of tem- 
perature with time is given by P] 



dT _ TV (T) - V (T) df / 87r3g^ 



1/2 



rp6 



(4) 



where Mpi = 1.22x lO^^Gey is the Plank mass. The first 
term is the contribution of reheating. It describes the 
increase of energy density of the plasma due to the release 
of latent heat during the phase transition. The second 
term is just — iJT; it accounts for the decrease of energy 
density caused by the expansion of the universe. Before 
and after the phase transition df /dt = 0, and Eq. (^) 
gives the well known relation t — £^Mpi/T'^, with ^ ~ 
0.03. During the phase transition, the coupled integro- 
differential equations (|^,^ must be solved numerically. 

To evaluate the effect of a changing wall velocity on 
baryogenesis we need to know how the baryon production 
depends on Vw As the bubble wall sweeps through space, 
it leaves behind a baryon density 
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-cr„ 



dx , 
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where F^.s ~ 20a^T |jT^ is the weak sphaleron rate, and 
ni is the net left handed density in front of the wall. 
The exponential accounts for sphaleron relaxation of the 



baryon asymmetry for small velocities. The coefficient c 
depends on the squark spectrum and is ~ 10. 

The left-handed density can be assumed to be of the 
form UL {x) = Ae""™^/^ , where D - 100/T is an ef- 
fective diffusion constant for the chiral asymmetry [ pO[ . 
The constant A depends on the CP violating force at the 
bubble wall that sources the asymmetry. It is in general 
proportional to the wall velocity, A oc Vw, so integration 
of Eq. (|) yields 



C 



+ cT^sD/v 



(6) 



where C does not depend on v^. This analytic approxi- 
mation describes qualitatively the dependence of ns on 
Vw It has a peak at Wpoak = y/cT^gD ~ 0.02. A cal- 
culation of the coefficient C is out of the scope of this 
paper, since we will only be interested in relative values 
TT-B (vw) /nB (vo)- Note that if the initial wall velocity is 
''^o 3> Wpcak, then a decrease of will produce an en- 
hancement of the BAU, since in that case ub ^ ''^w^ ■ 
the contrary, for small velocities '^-^ Wu,, so a velocity 
decrease will cause a suppression of the final BAU. 

As varies during the phase transition, ub {vw {t)) 
gives only a local baryon density, generated in a vol- 
ume dV (t) = — Vrotai/ {t) dt. The final baryon den- 
sity is the average over the expansion of bubbles |2^] 
B — ^ J Ub {t)dV. According to Eq. (^, this is re- 
lated to the result obtained with a constant wall velocity 
vq hy B — Sub (vq), where the factor 



df yp + ^'gcak/^O 

dt Vn, (t) + vl^^Jvw (t) 



dt 



(7) 



gives the enhancement or suppression due to the effect of 
reheating. 

We consider three sets of parameters for the elec- 
troweak phase transition. We take as case A the SM with 
an unrealistically small value of the Higgs mass, which al- 
lows for a sufficiently strong first-order phase transition. 
Therefore, we set the values D — 0.2 and E ~ 0.006 
0, and we choose A — 2E in order to fulfill the con- 
dition V (Tc) /Tc ^ 1, which is required for avoiding the 
washout of the BAU. For the friction of the plasma in 
the SM we assume the rough value 77 ~ 1 p-Q] . 

For the MSSM, the one-Higgs potential (|l[) can be used 
as an approximation in the case in which only one Higgs 
boson is light . In this scenario the parameter E can be 
at most an order of magnitude larger than in the SM. Al- 
though this limiting situation is hardly achieved in prac- 
tice, we can use it to simulate, without departing from 
reasonable values of the parameters, an extension of the 
SM that is most favorable for baryogenesis. We thus 
choose for case B the values E — 0.06 with A — 2E, 
and D ^ I. In this case we assume a friction coefficient 
77 ~ 70, in accordance with recent calculations for the 
MSSM with a light right-handed stop fll. 
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As we use it only to compute the dynamics of the phase 
transition, V {(j), T) need not be a real perturbative effec- 
tive potential. Instead, we can regard the field (p as an 
effective order parameter, and use Eq. (|l|) to model the 
phase transition dynamics |^,|f^ ■ In that way the param- 
eters D, E, and A can be chosen so that the free energy V 
has the same thermodynamical properties of the theory 
we wish to study. The relevant quantities are the latent 
heat L defined above, the surface tension of the bubble 
wall, cr = J {d(p / dx)'^ dx\Ta , ^^'^ the correlation length, 
given by = d'^V/d(j)'^\yi^T^-f T^. In our model these pa- 
rameters are given by L/T* = 8D {E/Xf (l - E'^/XD), 
ajTl ^ 2y/2E^/2,\^/^, and C^c = y/W^jE. We thus 
choose the parameters of case C in accordance with recent 
non-perturbative lattice simulations of the MSSM in the 
light right-handed stop scenario, which give i/T* ~ 0.4, 
ajTl ~ 0.01, and ^T^ ~ 5 ||. This case represents a 
more physical situation, as it allows electroweak baryoge- 
nesis for experimentally viable values of the Higgs mass. 

In Fig. 1 we have plotted the evolution of the phase 
transition for case A. We have defined a dimensionless 
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FIG. 1. Dimensionless temperature a, fraction of volume 
in the symmetric phase /, and bubble wall velocity , as 
functions of dimensionless time r, for the parameters of case 
A, with and without including the release of latent heat. The 
vertical lines delimit approximately the phase transition in- 
terval. 

temperature a = (T — Tq) / (Tc — Tq), and a dimension- 
less time T — (t — tc) I [Iq — tc), where t^ — S^Mpi/T^ is 
the time at which the universe reaches the temperature 
Tq if reheating is ignored. The result obtained when ne- 
glecting the reheating is plotted with dashed lines. In 
that case the duration of the phase transition is so small 
that all the parameters that enter the BAU can be ap- 
proximated by their values at the temperature T„ of the 



onset of nucleation ||15| . When the effect of latent heat is 
included, the plasma heats up after the beginning of nu- 
cleation because the expansion of the universe does not 
remove energy fast enough to compensate the release of 
latent heat. The bubble wall velocity thus decreases, and 
the phase transition slows down until an equilibrium tem- 
perature is reached, at which all the released latent heat 
goes into expanding the universe. Finally, the phase tran- 
sition completes and the temperature decreases again. 
Due to the reheating the wall velocity approaches Wpcak, 
so the baryon asymmetry results enhanced. Integration 
of Eq. (|^) gives S* ~ 3. As expected, the enhancement is 
not considerable since, as the latent heat is not very large 
in this case, Vyj does not spend a long time near '^pcak- 

In case B, in contrast, the latent heat is much larger. 
The plasma heats up very close to Tc and the temper- 
ature remains constant for a long period of time (see 
Fig. 2). However, the BAU turns out to be suppressed 
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FIG. 2. Same as Fig. 1, but for case B. 



in this case. Since the initial wall velocity is less in the 
MSSM due to the larger friction of the plasma, the de- 
crease of Vw occurs in the region on the left of the peak. 
Eq. (Q) gives a suppression factor S* ~ 0.15. Note that the 
slow growth period, in which ~ 2 x 10^^ ~ 10~^?;o, 
gives only a contribution ~ 10~^ to S. The baryon 
asymmetry is essentially generated in the initial stage 
in which bubbles fill a 20% of the volume of the uni- 
verse. There exists an additional suppression with re- 
spect to the maximum BAU, since the initial wall veloc- 
ity vq is a little below Wpeak- The total suppression factor 

is 2 {Vo/Vpcak + Upeak/'i'o)"^ S ~ 0.1. 

The results of case C are plotted in Fig. 3. Here the 
latent heat is smaller than in the previous case, so the 
final velocity is larger, v^, ~ 5 — 7 x 10"^, and the phase 
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transition occurs in a shorter time. As a consequence, the 
suppression of the baryon asymmetry due to reheating is 
less severe, S* ~ 0.4. 
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FIG. 3. Same as Fig 1, but for case C. 

Finally, it is important to notice that in some cases the 
simplifying approximation (^) may fall short of describing 
the behavior of (wu,), especially for small Vw For such 
cases, further suppression may arise in a more rigorous 
calculation. For example, in Ref. a crossing through 
zero of ub occurs at t;^, ~ 10^'^. That happens because 
the CP violating force changes sign near the bubble wall, 
and so does the chiral asymmetry. This gives an opposing 
contribution to n^, which becomes important for small 
wall velocities. Hence, in our last two cases baryon num- 
ber densities of opposite sign would be generated in dif- 
ferent regions of the universe. The negative contribution 
of the slow growth period to Eq. (0) can be roughly esti- 
mated by multiplying the ratio ub (10~^) /ub (10~^) by 
the fraction of space A/ that is filled during this stage. 
According to the results of Ref. , that ratio depends 
strongly on the bubble wall width and on the gaugino 
mass parameters. If we take for instance a value ~ —1/3, 
then in case B, with A/ ~ 0.8, this contribution is of the 
same order of the previously generated baryon asymme- 
try. In case C, with A/ ~ 0.4, the suppression factor 
would decrease a 30%. 

To summarize, we have calculated the decrease of the 
bubble wall velocity that occurs during the electroweak 
phase transition as a consequence of reheating, as well 
as its effect on baryogenesis. In the MSSM this effect 
tends to suppress the baryon production. For the light 
stop scenario we have found a suppression factor S ~ 0.4. 
Although this is not severe, we must stress that S can be 
smaller if the baryon density changes sign for small wall 
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